(N 



•3 



An inversion formula for the spherical mean 
transform with data on an ellipsoid in two and three 

dimensions 



^ ■ Yehonatan Salman 

Bar-Ilan University 

< 

Abstract 
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In the articles [1] and [2] of D. Finch, M. Haltmeier, S. Patch and D. Rakesh 
■ inversion formulas were found in any dimension n > 2 for recovering a smooth 

function with compact support in the unit ball from spherical means centered 
on the unit sphere. The aim of this article is to show that the methods used in 
[1] , [2] can be modified in order to get similar inversion formulas from spherical 
means centered on an ellipsoid in two and three dimensional spaces. 
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1 Introduction and mathematical background 

in . 

Let n be an integer equals to 2 or 3. Denote by S n = {x G M. n : |x| = 1} 
the unit sphere in IR n . For a closed smooth bounded domain Q denote by dQ its 
boundary and by C£°(Q) the class of infinitely differentiate functions denned on 
W l which vanish outside Q. For any n positive real numbers ai,...,a n > denote 
by E = E ait ___ tan the standard solid ellipsoid in R n with axes a±, ...,a n , that is E = 
{x G M. n : a{ 2 x\ + ... + a^x^ < 1}. Let cISe be the measure on dE corresponding 
to the measure that comes from integration on the unit sphere, that is, in the two 



dimensional case we have 



'dE 

and in the three dimensional case 



/7T 
f(ai cos6,a 2 sin 
-7T 



PIT piT 

f(p)dSE(p)= / / /(oi sin cos 6 1 , a 2 sin sin 6 1 , a 3 cos^) sin 

'dE Jo J-TT 

The spherical mean transform R takes any function / G C°°(IR n ) to 

(Rf)(x,r)= f f(x + rco)dS(u) 



n-l 
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where dS is the standard surface measure on S n 1 . 

The goal of this article is to recover the function / G C^°(E) from the spherical 
mean transform Rf restricted to the surface dE x [0, oo) where E = E au _ >an and 
cl \ , . . . , are positive real numbers. That is, our goal is to recover / from spherical 
means with arbitrary radius centered on dE. We will give inversion formulas that 
recover / in the two and three dimensional case. Inversion formulas for the general 
case of an arbitrary ellipsoid in two or three dimensional space can be easily obtained 
from inversion formulas for the standard ellipsoid by using translations and rotations. 
The methods presented in this article are based on and generalize the methods pre- 
sented in [1], [2]. 

The inversion problem for the spherical mean transform has been studied deeply 
recently due to its applications in various mathematical, physical and scientific fields 
such as thermo and photoacoustic tomography, radar and sonar imaging, approxima- 
tion theory and other areas ([6, 8, 11, 12, 14]). The inversion problem has been found 
out to be extremely important in thermoacoustic tomogrpahy (TAT for short). In 
TAT a short duration EM pulse (electro magnetic pulse) is sent through a biological 
object Q which causes a thermoacoustic response from Q resulting in a propagation 
wave G(x, t) which describes the amount of energy absorbed at location x and time 
t > 0. The transducers (the receivers that measure the pressure wave G(x,t)) are 
placed on the boundary dfl of the biological object and from this data one has to 
recover the initial pressure wave f(x) = G(x,0) on Q. Recovering / can be useful 
in finding cancerous cells in the biological object Q since they absorb several times 
more energy than ordinary cells ([6, 15]). The recovering of the initial pressure wave 
is strongly related to the inversion problem for the spherical mean transform since 
the amount of EM energy absorbed at location x and time t > is the average energy 
that was generated on the sphere with center x and radius t at the initial time which 
is just the spherical mean transform (Rf)(x,t) of the initial pressure wave. 

2 Main result and known results 

From now on a±, a n will be fixed positive real numbers. The inversion formulas 
presented in this article differ in the two and three dimensional case. Both formulas 
are modifications for the elliptic case of Theorem 1.1 in [1] for the two dimensional 
case and of Theorem 3 in [2] for the three dimensional case. In the two dimensional 
case we have the following theorem 

Theorem 2.1. Let E = E ai>a2 C M. 2 be the standard solid ellipsoid with axes 01,02 
and f G Cq°(E), then for every x G E, 

fid 2 1 d 2 \ 1 f f°° 
f(x) = + -p-^j — ^ r(Rf)( P , r) log \r 2 - \x - p\ 2 \ drdS E (p). 

(2.1) 
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In the three dimensional case we have the following theorem 

Theorem 2.2. Let E = E ai ^ a2 ^ a3 C M. 3 be the standard solid ellipsoid with axes 
a±, 02, 03 and f G Cq°(E), then for every x E E, 

fid 2 Id 2 1 d 2 \ I [ f 00 

f{x) = -{^ 1 + ^ 2 + 4^)^J BE J r2(Rf)ip ' r)6 l r2 - lx ~ P|2 I drdSM > 

(2-2) 

where S is the Dirac delta function. 

It should be noted that the integrals that appear in the theorems above are well 
defined since / G C™(E). 

Inversion formulas for the spherical mean transform with centers on a sphere have 
been derived recently in [1, 2, 7, 10, 16]. In [2] an inversion formula was found for odd 
dimensions and later in [1] a similar inversion formula was found for even dimensions. 
In [7] a uniform inversion formula was found for any dimension. 

Inversion formulas with spheres centered on an ellipsoid have been derived quite 
recently ([4, 5, 9, 13]). In [9] an inversion formula was found for the two dimensional 
case with spheres centered on the boundary of an arbitrary smooth domain, the 
inversion formula recovers / modulo an integral operator which vanishes for elliptical 
domains. In [5] a generalization of the results obtained in [9] was found for any 
dimension. A different inversion formula for elliptical domains for any dimension was 
found in [13]. The methods for recovering / presented in this article are different 
from the methods in [4, 5, 9, 13] and are based on the methods presented in [1] ,[2]. 
In fact the methods presented here are modified versions for elliptical domains of the 
methods in [1], [2]. 



3 Proofs 

We first prove Theorem 2.1. The proof of Theorem 2.1 relies on the following 
lemma which is a modification of Proposition 2.1 in [1] for the elliptic case. 

Lemma 3.1. Let E = E ai>a2 C M 2 be the standard solid ellipsoid with axes a±,a2- 
For every x,y G E,x 7^ y we have 

/ log I \x - p\ 2 — \y — p\ 2 \ dS E {p) = 27rlog \A(x - y)\ + C 

JdE 



where A 



ai 



a 2 

Proof. Start with the identity 



and C is a constant independent of x and y. 



\x — p\ 2 — \y — p\ 2 \ = 2 



x + y 



p) -(x-y) 
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which is true for every x, y,p G M 2 . Taking the log function on both sides of the last 
equation yields 



log||:r-p| 2 - \y-p\ 2 \ =log2 + log 

1 



x + y 



p) -ix-y) 



log 2 + log 



(M -\y\ )-P'(x-y) 



(3.1) 



Let us represent p in the form p = (ax cos 9, a 2 sin 9), then 

p-(x-y) = (ai cos a 2 sin #)-(a;i- 2:2-2/2) = ~yi) cos9 + a 2 (x 2 -y 2 ) sin9 

(cos 9, sin 0) • (ai(zi - 2/1), a 2 (^2 - 1/2)) = 



= (cos sin 1 



ai 



(xi - 2/1, rr 2 - 2/2) = e J • A(x - 2/) 



a 2 

Here • denotes the dot product. Hence from (3.1) and (3.2) we have 

1 



(3.2) 



log I a: -p\ —\y — p\\ = log 2 + log 



^(\x\ 2 - H 2 ) - e* ■ A(x - y) 



log 2 + log \A(x - y) \ + log 



\x\ 



M 2 _ e « 



A(x - y) 



2\A{x-y)\ \A{x-y)\ 
Integrating both sides of the last equation on dE with measure cISe(p) yields 



/ log \\x -p\ 2 - \y-p\ 2 \ dS E (p) 

JdE 

n | x |2 _ U|2 

27rlog2 + 27rlog|A(x-2/)| + / log 1 

J —IT 



A(x - y) 



2\A{x-y)\ 



\A{x-y)\ 



d9. 



Denoting e 1 ^ = A(x — y)/\A(x — y)\, using the change of variables 9 1— > 9 + ip and the 
2n periodicity of the function 9 >->■ e %e ■ e' 1 ^ = cos(9 — •?/>) gives 



/ \og\\x-p\ 2 - \y -p\ 2 \dS E (p) = 

JdE 



/7T 
log 
-7T 



X 



\v\ l 



Let us note that 



\x 



2\A(x-y)\ 
2 -\y\ 2 \A-\x + y) ■ A{x-y)\ 



cos 9 



d9. 



2\A{x-y)\ 



2\A{x-y)\ 



< l -\A-\x + y)\< 1 -\A- 1 x\ + 1 -\A- l y\ = 



l xj 4 1 yl vl 



< 1. 



Hence we can write 



\x\ 2 — \y\ 2 



cos a 



and thus 



/ 

JdE 



2\A{x-y)\ 
log | \x — p| 2 — \y — p\ 2 \ cISe(p) = 

log I cos a — cos 9 

-7T 

/7T 
log 

6» + a 



|d0 = 
9-a 



1 . + a 
- sin I sin , 

2 \ 2 V 2 



27rlog|A(x-y)| 



log 



sin 



/ lo, 



sin 



0-a 



d0. 



Using the change of variables t = 9 + a and t = 9 — a in the first and the second 
integrals respectively, and observing that the function t \ sin(t/2)| is 2tt periodic, 
we have 



/ log \\x -p\ 2 - \y -p\ 2 \ dS E (p) 

JdE 

log 

-7T 



dt. 



The important thing to note is that 2 f* log |sin (t/2)) \ dt is a constant C which does 
not depend on x or y, hence 



/ log \\x -p\ 2 - \y - p\ 2 \ dS E (p) 

JdE 

2nlog\A(x-y)\ + C. 



This proves the lemma. 



□ 



We now turn to the proof of Theorem 2.1 which is a modification of the proof of 
Theorem 1.1 in [1]. 

Proof of Theorem 2.1 : For / G C^>(E) and x e E let us compute 

/' poo 

/ / r (Rf) (p, r) log \r 2 - \x - p\ 2 1 drdS E (p) = 

JdE JO 

/ / /(y)log||y -p\ 2 - \x -p\ 2 \dydS E {p) = 
Joe Jm? 

/ f(y) \og\\y-p\ 2 - \x -p\ 2 \dS E {p)dy. 

Jr 2 JdE 
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Using Lemma 3.1 on the inner integral in the right hand side of the last equality gives 

r(Rf)(p, r) log \r 2 — \x — p\ 2 \drdSE(p) = 



dE JO 



J 

Jr 2 



f(y)(27r\og\A(x-y)\+C)dy. 



Hence by applying the operator a 1 2 d 2 /dx\ + a 2 2 d 2 /dx\ on both sides of the last 
equality we get 



id 2 id 2 

+ 



a\ dx\ ' a\ dx\ 



J / / r(Rf)(p, r) log \r 2 — \x — p\ 2 \drds{p) = 

J JdE JO 



Id 2 Id 2 ,. 

' f(y)(2nlog\A(x-y)\ + C)dy 



s a\ dx\ a?, dx\ 

271 L m (hjk + h 2 h) log |ai(xi ~ a2(X2 - m)ldy = 

2vr / f(y)(d 2 + d 2 )\og\a 1 (x 1 -y 1 ),a 2 (x 2 -y2)\dy = 
Jr 2 

2vr / f(y)2irS(ai(xi - yi),a 2 (x 2 - y 2 ))dy = 4n 2 f(x), 
Jm. 2 

where we have used the fact that (27r) _1 log \x\ is a fundamental solution of the Lapla- 
cian in M 2 . This proves formula (2.1) and thus Theorem 2.1. 

We now turn to the proof of Theorem 2.2. The proof of Theorem 2.2 relies on 
Lemma 3.2 which will be stated below. The method of the proof of Lemma 3.2 is 
taken from the first part of the proof of Theorem 3 in [2] . 

Lemma 3.2. Let E = E aiiCl2ia3 C M 3 be the standard solid ellipsoid with axes ai,a 2 ,as. 
For every x,y G E, x ^ y we have 



I 5 \\x -p\ 2 - \y -p\ 2 \ dS E (p) 

JdE 



7Y 



\A(x-y)\ 



a 1 
where A = \ a 2 
a 3 

Proof. By the coarea formula, if M is a two dimensional surface in M 3 , given by 
4>(p) = 0, with V0(p) 7^ at every point of M, then for every h G 



/ h(p)dS p = [ h(p)\V<l>(p)\5(<l>(p))dp (3.3) 

where dS p is a surface element of M. In [3] Chapter V it was proved that if dE is 
given by (f>(p) = a^ 2 p\ + a 2 2 p\ + a^ 2 p\ — 1 = 0, dS p is a surface element of dE and 
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dS E is the measure on dE coming from the unit sphere as defined above, then the 
following relation holds 

Using this and (3.3) with M = dE yields 

/ h(p)dS E (p) = — ^— / T^y { dS p = — ^— [ h(p)5(<f>(p))dp. 
JdE 010203 J d E |W(p)| aia 2 a 3 i M 3 

Hence in particular for h(p) = 5 \\x — p\ 2 — \y — p\ 2 \ we have 

f 6\\x-p\ 2 -\y- P \ 2 \dS E (p) = 

JdE 

/ S | \x - P \ 2 - \y - p\ 2 \ 5(a^ 2 pj + cq 2 p\ + ag 2 p\ - l)dp = 

/ $\\ x \ 2 ~ \v\ 2 - 2 P • (x - y)\ 5(a^ 2 p\ + a^ 2 p\ + a 3 2 p\ - l)dp. 
The change of variables Ai = a^pi, A 2 = ^2 V2, A3 = a^ps gives 
f 5\\x-p\ 2 -\y-p\ 2 \dS E (p) = 2 [ 6\\x\ 2 -\y\ 2 -2\.(A(x-y))\5(\\\ 2 -l)d\ 

JdE </R 3 

I ai \ 

where A = (Ai,A 2 ,A 3 ) and A = a 2 . Let Q be an orthogonal matrix 

V a 3 J 

satisfying QA(x — y) = \A(x — y)\e 3 where e 3 = (0, 0, 1). Using the change of variables 
A = Q~ l p gives 

f 6\\x-p\ 2 -\y-p\ 2 \dS E (p)=2 [ 5\\x\ 2 -\y\ 2 -2p-(QA(x-y))\5(\p\ 2 -l)dp = 
JdE Jm. 3 

2/ 5\\x\ 2 -\y\ 2 -2\A(x-y)\p,\5(\p\ 2 -l)dp. (3.4) 

Using (3.3) on equation (3.4) with the surface H given by <f)(p) = \x\ 2 — \y\ 2 — 2\A(x — 
y)\P3 — an d noting that |V0(p)| = 2|A(x — y) \ gives 

f S \\x-p\ 2 -\y-p\ 2 \dS E (p)= 1 / <5(M 2 -l)dtf. 

using the following parametrization for H 

\x\ 2 -\y\ 2 

H : ( Pl , P2 ,a), Pl ,p 2 G K, a = — rr, 

2|A(x-j/)| 

gives 

/ 5|k-p| 2 - |y-p| 2 |^ B (p) = — — ^ -/ 5(p 2 +p" + |a| 2 -l)dp 1 dp 2 . (3.5) 

i^h^ ~~ yji ^m 2 



So far we used (3.3) only for two dimensional surfaces but in fact formula (3.3) is true 
in any dimension and in particular it is holds for one dimensional surfaces. Hence 
using (3.3) on equation (3.5) for the circle C : 4>(p) = p\ + p| + \a\ 2 — 1 = and 
observing that |V0(p)| = a/ Ap\ + Ap 2 = 2y/l — \a\ 2 gives 



ft ||< - if - \v - v?\ ds E ( P ) = / -£ 



dC 7T 



i^/T^W \ A (x-y)Y 

This proves Lemma 3.2. □ 

Remark 3.3. As in the proof of Lemma 3.1 it can be shown that \a\ < 1 and thus the 
square root a/1 — \a\ 2 is defined well. 

We now turn to the proof of Theorem 2.2 which is a modification of the first part 
of the proof of Theorem 3 in [2] . 

Proof of Theorem 2.2 : For / G C yo (E) and x e E let us compute 



dE JO 



r (Rf)(p,r)5\r — \x — p\ \drdSE(p) = 
/ / f(y)8 \\y ~P\ 2 ~ \x -p\ 2 \ dydS E (p) = 

JdE JR 3 

/ f(y) 6\\y-p\ 2 -\x-p\ 2 \dS E (p)dy. 

JR 3 JdE 

Using Lemma 3.2 on the inner integral in the right hand side of the last equality gives 
/ / r 2 (Rf)(p,r)5\r 2 - \x - p\ 2 \drdS E (p) = / f{y) T j r ^- — -dy. 

JdE Jo Jr3 \ A {x-y)\ 

Hence by applying the operator a{ 2 d 2 jdx\ + a^ 2 d 2 jdx\ + a^ 3 d 2 /dx\ on both sides 
of the last equality we get 



Id 2 Id 2 Id 2 



a\ dx\ a\ dx\ a\ dx\ ) J dE J 



r 2 (Rf)(p, r)5 \r 2 — \x — p\ 2 \ drdS E (p) 



Id 2 Id 2 1 5 2 \ / w , 7T 



a 2 dx\ a\ dx\ a\ dx\ ) J^ \A(x - y) \ 

1 d 2 1 d 2 1 d 2 \ f £ . . 7T 

+ 77.7772 + -2 77-2 / f(y)TT7Z 77—77 7777777 77T\ d V 



a\dx\ a\dx\ a\dx\) J R3 \a 1 (x 1 - y ± ), a 2 {x 2 - y 2 ), a 3 {x 3 - y 3 )\ 
f(y)(d 2 + d 2 2 + dt)— -— -dy 



ai{x! - yi),a 2 (x 2 - y2),a 3 (x 3 - y 3 )\ 

-4tt 2 / f(y)5(a 1 (x 1 - y{), a 2 {x 2 - y 2 ), a 3 (x 3 - y 3 ))dy = -4n 2 f(x), 
Jr 3 

where we have used that fact that (— 4n\x — y|) _1 is a fundamental solution for the 
Laplacian in IR 3 . This proves formula (2.2) and thus Theorem 2.2. 
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